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Abstract 

The seesaw model of quark masses is studied systematically, focusing on its develop- 
ments. A framework allowing the top quark mass to be of the order of the electroweak 
symmetry breaking scale, while the remaining light quarks have much smaller masses, due 
to the seesaw mechanism, is presented. The violation of the GIM mechanism is shown 
to be small and the tree level FCNC are suppressed naturally. In this model, there are 
many particles which could contribute to the FCNC in the one-loop level. Parameters of 
the model are constrained by using the experimental data on K° — K° mixing and ex- 
The rare K meson decays K l g — > ifivv and K + — ► ir + ui> are also investigated in the 
model. In these processes the scalar operators {sd){v T v T ), which are derived from box 
diagrams in the model, play an important role due to an enhancement factor Mk/tti s in 
the matrix element < ir\sd\K >. It is emphasized that the Kl decay process through the 
scalar operator is not the CP violating mode, so B(Kl — > 7r°z/p) remains non-zero even 
in the CP conserved limit. The pion energy spectra for these processes are predicted. 



1 Introduction 

The idea of the seesaw mechanism for the neutrino mass [Q || was extended to quark masses 
sometime ago At that time, the top quark mass was not known and it was not expected 

to be as heavy as the electroweak symmetry breaking scale. Therefore all charged fermions 
were assumed to be lighter than 10 2 GeV. The motivation for the construction of the model 
was to give an explanation for the smallness of the masses of the charged fermions (quarks and 
charged leptons) as compared to the electroweak breaking scale. However, now we know that 



the top quark is as heavy as O(10 2 GeV) and that has to be taken into account. Koide and 
Fusaoka (KF) || and T. Morozumi, T. Satou, M. N. Rebelo and M. Tanimoto |J studied the 
top quark mass problem in a seesaw model. Based on this successful approach, we can discuss 
implications of the seesaw model focusing on the experimental data such as flavor changing 
neutral currents (FCNC) and CP violation. In this paper, we present in detail the formalism 
and make a systematic phenomenological analysis of the implications of the model. A brief 
introduction to the Dirac seesaw scheme is given in section 2, in which the incorporation of 
the top quark is outlined. In section 3 we show how the mass hierarchy and the flavor mixing 
are introduced. Section 4 is devoted to estimating tree level FCNC. In section 5, the flavor 
mixing in the charged currents is studied and the parameterization of the generalized Cabbibo 
Kobayashi Maskawa (CKM) matrix is specified. [[UJ Jl7j In section 6, the K° — K° mixing and 
the CP violating parameter ex are studied and the constraint on the parameters of the present 
model is obtained. In section 7, the rare K meson decays K — > i^vv are studied in the present 
model. The effect of the scalar interaction on the resulting pion energy spectrum is given. Our 
conclusions are presented in section 8. 



2 Seesaw Model and Top Quark Mass 

Let us first summarize the ideas of the seesaw model for quark masses. The gauge group of the 
standard model (SM) is extended to SU{2)l x SU{2)r x £7(1), so that the usual Yukawa mass 
terms for the charged fermions are not allowed, since now the right handed charged fermions 
belong to an SU{2) R doublet. That is: 



<W(2) fl {^lHr) ^ 0. (1) 

The SU{2) R gauge boson must be much heavier than the SU{2) L gauge boson because we 
do not see any deviation from the V-A structure of the charged currents, thus the symmetry 
has to be broken at a sufficiently high energy scale. The representation of the Higgs breaking 
SU(2)r is chosen in such a way that no renormalizable mass term for the standard-like fermions 
is allowed. Therefore, for example, the Higgs bi-doublet M(2, 2) under SU{2)l x SU{2)r is 
excluded. A single Higgs doublet, 4>r(1, 2), is used for breaking SU(2)r. With this Higgs field, 
the possible mass term for the standard-like quarks is the dimension 5 operator, 

1 -IpLtLfyR, (2) 



A-NEW 

where Anew is some new physics scale. The scale must be larger than rjn (VEV of 4>r). Quarks 
and charged leptons acquire their masses, which are much smaller than the electroweak breaking 
scale, t]l (i.e., VEV of 4>l)i through the dimension 5 operator: 

m f = 7 ' ( 3 ) 

^NEW 

where -r- 25 — << 1. Now we can think of what is the fundamental theory behind this dimension 

A-NEW 

5 operator. The interaction can be reproduced if we introduce new heavy isosinglet quarks (and 
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charged leptons) and integrate them out. These isosinglet quarks can have bare mass terms, 
and the new physics scale, Anew, is identified with that mass. Such a renormalizable theory 
is given, in the quark sector, by the following Lagrangian terms: 

C = -y L ipL<f>LU R - y R ?p R <p R U L + (h.c.) - M V UU. (4) 

Using this Lagrangian, we may compute the Feynman diagram in which the heavy singlet quark 
is exchanged, obtaining for the dimension 5 operator: 

1 _VLlf R (5) 



Anew Mt 



However this formula leads to the suppression of the quark masses by a factor of SU(2) R 
breaking scale divided by Singlet quark mass compared to the electroweak breaking scale, hence 
it cannot be applied to the top quark. In order to prevent the seesaw mechanism to act for 
the top quark we assume that the corresponding singlet quark has a bare mass much smaller 
than the SU(2) R breaking scale. In this case, ignoring flavor mixing for the moment, the top 
quark only acquires mass through a dimension 4 Yukawa coupling, instead of the mechanism 
of exchange of a heavy singlet quark. To illustrate, let us consider the extreme limit in which 
the bare mass is set to zero in Eq.(|J), so that we have 

£ = -VlT)i£lTr - y* R riRT L t R + (h.c). (6) 

These are dimension 4 operators and the singlet quark cannot be integrated out because the 
bare mass term is absent. The diagonalization is performed through maximal mixing for the 
right-handed sector, as follows: 





(7) 

R \ " / R 

(the superscript m denoting the mass states), leading to: 

m t = \y L \VL, m T = \y R \r]R, (8) 

and mt is of the order of the experimental value. If we retain the bare mass term, i.e., —MtTT, 
the mass formulae are given by 

i I. i VlVr 



m T = J\y R \ 2 v 2 R + M T . (9) 
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3 Quark Mass Hierarchy and Flavor Mixing 



We will now extend our analysis to the case where flavor mixing is present and the number of 
generations is three, based on the following Lagrangian: 



-y l i LE L^ L E R - y^ RE L^cj> R E{ + h.c. 
-WMiU* - D i Mi D D i - E i MgE i 7 



(10) 



where My, Mp and M| are real parameters. The representations under the gauge groups of 
ordinary quarks and charged leptons are assigned as: 



Uf) 



(2,1,1/3), 



rOi 
L 



Ur 



;i,2,l/3), L 



Qi 



U 



R 

L,R 




: (2,1, 



: (1,2,-1), 



(1,1,4/3), Di jR :(l, 1,-2/3), E 



L,R 



1,1,-2 



We do not introduce singlet neutrinos and so there are no tree level neutrino masses. Let us first 
study the charged lepton sector and the down quark sector. The masses of the charged leptons 
and down quarks are much smaller than O(100GeV), so they can be treated in a unified way. 
We consider first the down-like quarks. In order to find the mixing angles (CKM matrix) among 
both left and right chiralities, it is convenient to perform unitary transformations among the 
ordinary quark fields such that the singlet-doublet Yukawa couplings, i.e. y L and y 0R , become 
triangular matrices 0, 

/ 2/i 



(12) 



2/2i 2/2 

\ 2/31 2/32 2/3 

where U is a unitary matrix. The diagonal elements in Eq. ([l2D are real and the off-diagonal 
entries are complex. We perform this same type of transformation in the charged lepton sector. 
Accordingly, we obtain new bases for the SU(2) doublet quarks and leptons, 



< 


= U j UL u 0L , 


u'r 


= U RU u 0R , 


d' L 


= U DL d L, 


d' R 


= U DR d R, 


e' L 


— U EL e 0L , 


e' R 


— U ER e 0R , 


i 


= U EL v 0L , 


v 'r 


= U ER V R. 



(13) 
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We note that the transformation on neutrinos is chosen to be the same as that on charged 
leptons. As a result, the down quark and charged lepton mass matrices are written as follows: 



M 



v 



Mi 












yLDlVL 

















yhD2\r\L 


yLD2VL 














yLD3lVL 


yLD32VL 


yLD3VL 


VrdiVr 


V*RD2lVR 


y rosier 


M D 











VRD2VR 


V*RD32VR 





M s 











yRD3VR 








M B 











yiEiVL 

















yLE2lVL 


yLE2VL 














yLE3lVL 


yLE32VL 


VLE3VL 


VreiVr 


V*RE2lVR 


y*RE3lVR 


M E 











URE2VR 


y*RE32VR 





Mn 











yRE3VR 








M T 



(14) 



(15) 



These mass matrices can be approximately diagonalized by using the following unitary matrices, 

VqDL = 



Vn 



ODR 



1 yLDVL/M v 

■l/M v y LD ^ VL 1 

1 yRDVR/Mv 

-l/M v y RD ^ R 1 



(16) 



where M 



diag(Mr>, Ms, Mb). As a result, Mt> is approximately diagonalized as 
/ - 



VqdlMvVqdr 



\ 



-yLDi^yRDi 
-yLD2i n i^-yRDi 

-yLDSl^VRDl 







-yLDlV^y R D21 -ED1^&31 



-yLD2 n §fyRD2 



-yLD2 n §fy RD3 2 



-VLD32 m $fVRD2 -yLD3 ll §fyRD3 



o 













































M D 













M s 













M B 


) 



(17) 

Approximate eigenvalues of the down quark masses are given by the diagonal elements of 
Eq.flT7|) as: 



For charged leptons, similar results are obtained. 



m b 


r\J 


Vr 

yLD3yRD3VL, 

m b 


m s 




Vr 

-T7-yRD2yLD2f]U 
Ms 


m d 




Vr 

w yRDiyLDlVL, 

Md 


m D 




M D , 


m s 


r\J 


M s , 


m B 




M B . 


s are 


obtained. 



(18) 
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3.1 Top quark mass and hierarchy of up type quark masses 



Now we turn to the up type quark masses. In the previous section we saw that the top quark 
mass, in the absence of flavor mixing, is given by Eq.(|9|), rather than the seesaw type formula. 
The extension to the case where flavor mixing is present, and Mt is not exactly zero, but 
restricted to Mt « T]r, was done in |§. Here we extend this previous analysis in a way that 
can be applied to the case where Mt is as large as t]r. First we derive the formulae for the top 
quark mass as well as the other light quark masses by solving the eigenvalue equation. Consider 
the up mass matrix: 

Mu={ ° t VLV J LU ), (19) 

where Mu = diag(Mu, Mc, Mt)- The corresponding eigenvalue equation for the quark masses 
is given by 



det(MuM j u - A) = det 



VlyOLUVoLU ~ A 



V0LuM U r] L 



VLMuyl LU VrvIruVoru + M u - A 



0. 



(20) 



The equation which determines the eigenvalues of the order of rji 2 (or smaller than rji 2 ) is 
reduced to a cubic equation: 



det 



Volu {l - Mu^yomj^VoRU + M u) 1 M U } yl LU - A 



0. 



(21) 



where we use the normalized eigenvalue A = -y. We further use the following expansion: 



[—AM 2 —) + (—AM 2 —) 2 
V M Mq v M ilV 



— , (22) 
M ' 1 ; 



Mn 



Mr 



M, 



c 



vUVoRuVoRuhs + Mr 



(23) 



Note that in the traditional treatment of the seesaw model, Eq.flzz]) is expanded by the inverse 
power of the singlet quark mass matrix, i.e., diag(Mu, Mc, Mt)- However the smallness of Mt 
compared with t]r does not allow us to do so. M is chosen in such a way that the expansion 
by the inverse power of M is still regular even in the limit of vanishing Mt- By keeping the 
dominant coefficients of X n (n = 0, 1, 2, 3), the eigenvalue equation becomes, 



F(A) 



= -A 3 + {c 33 n 33 x T 2 }\ 2 

— {£33^33 {£ 22 ~ ^23^32/^33) (72-22 ~ 72 2 3723 2 /7233)X T ' 2 X C ' 2 }A 

+ X T 2 X c 2 X u 2 det£det1l = 0, 



(24) 
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where C tJ = {yl LU y 0L u) ip = {ylrnVmu)^ *u = j^,X c = jjfc, and X T = -j=Ja==^ . By 
noting that X\j << Xq « Xt = 0(1), the approximate solution of Eq. (|24D is 

VlVr 



m t = \/C 33 n 33 - 

Im t 2 + n 33 r) 2 R 

m c = J(C 22 - C 23 C 32 /C 33 ) (TZ 22 - Tl 23 ll 32 /Tl 33 ){'^j^-} , 

Mc 



m = / detC detU rWtei (2B) 

V (£22 - C 23 C 32 /C 33 )(n 22 - K 23 K 32 /K 33 )C 33 K 33 X M v 1 ; 

The heavier eigenvalues are also obtained as follows: 

mu = Mu, 
m c = M c , 



m T = \JM T Z + TL 33m 2 . (26) 

To compare Eq.fl2"5p and Eq.fl2"ED with the formulae obtained in the absence of flavor mixing, 
Eq.(^), we need to explain a little bit more. The mass formulae for the light quarks are obtained 
in a weak basis in which the singlet quark mass matrix M is diagonal, and the singlet-doublet 
Yukawa couplings are general. However, the formulae are simplified if we go to the weak basis 
in which the singlet-doublet Yukawa couplings are in the triangular form. For example, 



£33 — VLU3, V "^33 — VRU3 



C 22 — C 23 C 32 1 '£33 — yLU2, Y^22 — 7^23^32 — VRU2, 

VdetC = (yLui)(ym2)(vim), VdetK= {ymx) (yjora) {vruz) • (27) 

In this basis, the light quark masses are determined by the diagonal elements of the singlet- 
doublet Yukawa coupling: 

VlVr 



m t = {yLm){yRm) 



Mt 2 + y 2 RU3 V R 



, w srVLVR, 

m c = {yLU2){yRU2){-ri—}, 

m u = (yLui)(Vmi){j^}- (28) 

The top quark mass in Eq. (p8|) reduces to the result obtained in Eq.@. Therefore, we conclude 
that in the triangular basis, the off-diagonal elements of yu and y R can be safely neglected, since 
ignoring them does not affect the mass eigenvalues significantly. In the following section, it is 
shown that they are related to the size of the FCNC. Finally, we check the formulae Eq.(p6|) 
and Eq.(p8|) by showing that detM. u = m u m c m t mumcmT ■ In fact, detM.u is given by, 



detMu = si 9n(p)M U iP(i)Mu 2 p( 2 )M U3 p( 3 )Mu^)Mu5P(5)-M 



U6P(e) 



p 



MuiAM U25 MumM U AiM m2 M m3 

(yLui)(yLU2)(yLU3) {ym\){VBm){vwa) (VlVr) 6 - (29) 
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With Eq. (|26|) , Eq.(f28|) and vtlt = \J Mt 2 + UrusVr, detAiu is equal to the products of mass 
eigenvalues of the six quarks. 

4 Tree Level Z FCNC 

In this section, we study the tree level FCNC due to neutral gauge boson exchange. We first 
derive the theoretical form of these flavor changing currents. This will show how they are 
naturally suppressed, and enable us to perform quantitative calculations. Using the estimates 
and the present experimental bounds on rare B and K decays, we examine whether the tree 
level FCNC can significantly contribute to them or not. We can also compare the tree level 
FCNC effects with one loop GIM suppressed contributions of the standard model. 
Let us start with the relevant part of the Lagrangian: 

M 2 Z 2 M 2 2 Z 2 2 , N 

C = Z 1 J 1 + Z 2 J 2 + + (30) 

where Mi and M 2 are the mass eigenvalues of the neutral gauge boson mass matrix, 

/ cos 2 9 W n \ / sin 2 8 W tan 2 9 W tan 2 6 W \ 

M 2 — M 2 I cos26W U I _i_ M 2 I cos26» w - v / cos2^7 1 /q-i\ 

m - ivi WR \ 0' WL I tan 2 e w 1 I ' V 6L ) 



Vcos 28yy tan 2 9w 



where Mwr and Mwl are masses of the charged SU(2) gauge bosons and singly = e/g with 
SU(2) gauge coupling constant g and U(l) em gauge coupling constant e. 

/ cos 9 sin 9 \ „ r n I cos 9 — sin e \ ( m 2 2 o \ , ooN 



sin# cos 9 / I sin# cos 9 \ M 



When the breaking scale of SU(2)r is much larger than that of SU(2)l, i.e. Mwr >> Mwl, 
the eigenvalues and the mixing angle are approximately given by: 

m2 _ Mwl 2 M WL 4 tan 4 9 w 



M 2 



cos 2 9 W M W r cos 2 9 w 
2 M W r 2 cos 2 9 w 



cos 29w 



2 



2 Vcos 29 W 2 n M WL 
tan 29 = — tan 6W 5 , (66) 

cos 2 9 w M W r 



and Ji and J 2 are written as: 

Ji = Ji°cosfl- J 2 °sin#, 
J 2 = J 2 °cos9 + Vsinfl, 

Ji° = - 1 t^(-J3L + sin 2 9 w Jq), 



cos 9 



J 2° = g-r==t={-Jm - tan 2 9 w Jzl + tan 2 9 w Jq), (34) 
Vcos 26% 
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where J^l, J3R, and Jq are the SU(2)l isospin current, the SU(2)r isospin current, and the 
electromagnetic current, respectively. They are defined as: 



Jsl, = E Jsr, = E i>R 0i i^i>R 0i , Jq, = E ^ W 0a (35) 

i=l "' i=l ' a=l 

We write the isospin currents in terms of the mass eigenstates, 

= E V&dZ, 4* = E V™ R d R , (36) 

a=l a=l 

where Vdl an d Vdr are the 6x6 unitary matrices that diagonalize the seesaw type mass matrix, 

T/ f I 1/0LD% \ T/ _ ( m d \ 

^ i Wite M C J ^ - [ m D J • 
Now the neutral currents in the down quark sector are written as: 

■hh = --^22 DLaf3 d L a 'j fJi d L f3 , 

1 a/3 

J3R = - 9E 2 '-D Rapd R° '7 'pd R^ (37) 

where 

3 6 

- 2^ ^DL *DL ~ d _ 2^ DL V DL, 
i=l 7=4 

Z>DRa.(3 = j:V&V* a = S*-j:v£V&. (38) 
i=l 7=4 

In this way the tree level exchange of neutral gauge bosons gives rise to the effective AF = 1 
Lagrangian, 

C AF=1 = V2G F E {ZDLc0d?'Y l *d L ff rt>fv L i + f3Z DRa ^ lii d R ^^v R % } 
V2G F E {ZoLapd^l^ih^h 1 - 2 sin 2 6 w T^r) 

+ (3Z DRa ^ lfl d/(T^YlR - 2 sin 2 OwF-fl')}, (39) 

where the coefficient (3 is defined as (3 = M^/M^ . As we discuss below, the strength of the 
tree level FCNC couplings Z DR are enhanced by a factor of 1/(3 compared to Z DL . Therefore 
the AF = 1 FCNC due to Z x is of the same order of magnitude as that of Z 2 . 
We now turn to theoretical estimate. We first show how the FCNC among the ordinary down 
type quarks is suppressed. It turns out that the suppression factor is m£mj jMw R for Zdl 
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and mdmJ/Mw L 2 for Zdr- As seen in Eg . (|39D , the FCNC couplings among the ordinary 
quarks are given by: 



ZoLij — —V D % [vJ L , 

Z DRij = -VfXt (40) 



In order to estimate the FCNC using Eq. flioD , we must know Vdl and Vdr- F° r this purpose 
we can start with an approximate parameterization obtained from the diagonalization of the 
seesaw mass matrix and unitarity (see Eq . fll~2|) and Eq. fll3|) ), then we get 

VDL ~ { 1 J VoDL ~ { -l/M v y 0L J VL U DL 1 ) ' 

VDR - { 1 J W - ^ -1/M C W^ D « 1 J ' (41) 

where Udl and ^7d.r are unitary matrices that approximately diagonalize the 3x3 effective 
light quark mass matrix: 

- Udl (yoLD 1 j^yi RD ) U DR = m d . (42) 



v 



Notice that with a suitable choice of the unitary matrix U, an arbitrary matrix yo can be 
transformed into a triangular matrix, as in Eq. (|T2"D . Suppose U DL y Q DL and U R y oR are such 
triangular matrices. We immediately find that the effective light quark mass matrix is approx- 
imately diagonalized, 



TT \ I VLVR f \TT 

U DL [VOLD ~J^yORD ) U DR 



( VDLIJJ^VDRI yDLljhy*DR21 yDLlJJ^y* D R31 \ 

yDL21JJ^yDRl yDL2J^yDR2 yDL2Jj^y*DR32 
\ yDL31JJ^yDRl yDL32J^yDR2 yDLZ-j^yDRZ J 



VlVr, (43) 



where we assume the following hierarchy of the singlet quark masses, 

M D > M s > M B . (44) 
Therefore the light quark masses are approximately given by the diagonal elements of Eq.( 



m d = ~ [ yDLi^—yDRi ) VlVr- (45) 



Vi 



Substituting the approximate parameterization of Vdl and Vdr into Eq. fl4T|) , we obtain the 
following formulae for the FCNC: 



Z DL = -{U D LyOLDY^2yi LD UDL}ij 

= -m d l { U DR y mD y ] mD UdrY 1 ^ 3 /r] R 2 , 
z dr = -™d{U DL yvLDy\ L DUDLY X id mJ h L 2 . (46) 
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We must check that the right hand side does not necessarily vanish when % ^ j. We can see 
that this is indeed the case because UlUol and U R yo R are triangular matrices, and their off 
diagonal elements, y^x and y^ R (i ^ j), are non-zero in general. To leading order the off 
diagonal elements are given by: 

(i<j), 



yRDjy 2 RDt 

V *DLj 



(i<3), 



{u ] D mL D yl LD u DL }^ ~ { ) J !; (47) 



so the FCNC couplings are suppressed as: 

Z D Lij = -{m d m d 3 /2M WR 2 }g 2 (U ] DR y RD y RD U DR )^ 1 

yRDi \yHDj J 

Z DLij = Z* DLji (i> j), 

Z DRij = -{m d l m d i/2M WL 2 }g 2 (Ul )L y LD yi D U DL )r j 1 
~ (m d W/2M WL 2 )^^ f-O (i<j), 

yLDi \ yiDj J 

Zona = Z* DRji (i>j). (48) 

Therefore the tree level FCNC among the light quarks is naturally suppressed by a factor of 
(quark masses) 2 divided by an (577(2) breaking scale) 2 . Also, they are proportional to the off- 
diagonal elements of y^-. If all of these vanish, there are no FCNC among light quarks. 
We can now confront our result for the FCNC with the present experimental bounds from rare 
K and B decays. 



'Rsd 



'Rbd 



m s m d 


yhsd 


2M WL 2 


yid 


m b m d 


yhbd 


2M WL 2 


yhd 


m b m s 


yLbs 
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>VLs, 



'Rbs — „, 2 



3.1 x 1(T 7 {^ 

2y L d \yLsJ 

±) 2 = 7.8 x io-6{^ (^-)\ 
VLb) 2y Ld \y Lb J 

9 ^ - 1.6 x 10" 4 {^ ( ±)\ 



'Lsd — 



>Lbd 



■'Lbs 



2M WL y Ls \y L b) lyhs \yhb, 

g m s m d y Rsd ( _9_Y = 12x iq- 8 ( 400 V rVRsd ( 
P 2M WR 2 yRd \y Rs ) ~ • \M WR (GeV)J X 2y Rd \y Rs ) h 

m b m d y Rbd ( _g_\ 2 = 31x 1Q - 7 ( 400 \\y Rbd ( _9_\\ 
P 2M WR 2 y Rd \y Rb ) " ' ^ \M WR (GeV) J X 2y Rd \y Rb ) h 

m b m s y^ fj_\ 2 ( 400 \ \vn»(_g\\ 

P 2M WR 2 yRs {y Rb ) \M WR {GeV)J X 2y Rs \y Rb J h 
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p 



(49) 



where we use m d = 10 MeV, m s = 0.2 GeV, m h = 5 GeV and M Wl = 80 GeV. We also use the 
following notations. Zu, s = ZDL32,VRbs = VRD32 and so on. The tree level contribution to the 
branching ratio B — * X s e + e~ and K + — > n + vV is given by, 



(50) 
(51) 



B(K + - 


> 7r+z/z/)| Tree 3 




B{K+ 


-> 7r°e+z/) 2 


W I 2 



8(5 -> X,e+e- 



|Tree 



X c /-77) 



P,4 U 



2 sin 2 ^) 2 + 4sin 4 ^} 





2 


+ /3 2 


ZRsb 


2 




V cb 


2 



where P c = 0.538 and P s = 0.986 are the phase space factors for the charm quark and the strange 
quark, respectively. Using B(B -> XjTV) ~ 10%, |V C £ KM | = 0.04, -> 7r°e+z/) ~ 5% and 

| V™| = 0.22 [|TJ, the order of the magnitude of the tree level contribution to these processes 



arc, 



B[B - 
B(K^ 



X s e + e~)\^ ce < 1.4 x KT 9 , 

^-16 



n + uV)\ Tiee < 4.5xHT iD , (52) 

where M WR > 400 (GeV) and the combination of the coefficient of Yukawa coupling and gauge 
coupling in Eq.(|49|) is set to be 0(1). The experimental bound |U| and the prediction of the 



(53) 



standard model |14| of B — > X s e + e are, 



B{B ->X s e + e-)| Exp . < 5.7 x 10" 5 , 
B(B -> X s e + e-)| SM ^ 8.4 ± 2.3 x 10" 6 , 



while the recent measurement of K + 



7r + isv 



I3| and the prediction of the standard model 



18] are given by, 



B{K + 
B(K^ 



7T + zV)|e XP . 
» 7T + Z/77)| SM 



4.2^ x 10- 10 , 



0.6 ~ 1.5 x 10 



-10 



(54) 



Compared with values in Eqs.(|52D, fl53|) and (|54|) , the tree level FCNC is found to be negligibly 
small. In section 6 and 7, we investigate the FCNC in one-loop level. 



5 The flavor mixing in the charged currents 

In this section, we present the approximate parameterization of the flavor mixing in the charged 
currents which is analogous to CKM matrix in the standard model (SM). The difference between 
the CKM in the standard model and that of the present model is that the flavor mixing in the 
present model is 6 x 6 rather than 3 x 3 in the SM. There is also right-handed flavor mixing in 
the present model. As we show in the following, the flavor mixing consists of the singlet-doublet 
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mixing and the flavor mixing among doublet quarks. The charged currents of the present model 
are written as, 



Jl» = UL^V^dL 13 , 

where V L and V R are the generalized CKM matrix and are given by 



(55) 



V, 



H V ULai V DLi/3 ~ V 0ULai U tj V 0DLj/3, 
i i 
3 3 

H V URai V DRi/3 = Y V 0URajUij V 0DRj/3, 
i i 

where Ufj and Uff are 3 by 3 unitary matrices and are defined by, 



a/3 



V R 

V aj3 



(56) 



V) 



U L = U UL U DL , 

u R = u UR u DR . 



(57) 



Uul, Ujjr, Udl, and Udr are 3x3 unitary matrices which transform the singlet-doublet Yukawa 
matrices youL,youR,UoDL, and y 0DR into the triangular matrices, yuL,yuR,VDL, and y DR . 
Voul/r(Vqdl/r) are matrices which diagonalize M.u and M.x> as, 



VoUL-MuVoUR z 
VoDL-M-vVoDR 




(58) 



where 





VurVr 


VdrVr 



yULVL 
yDLVL 



(59) 



Because the mass eigenvalues are not affected significantly by the presence of the off-diagonal 
element of the triangular matrices yuL,yuR, Vdl and yDR, it is legitimate to neglect the off- 
diagonal element of the triangular matrices. In the approximation, M. v and M. u are approxi- 
mated as, 



M u - 














ywiVL 






















yLmVL 






















ywsVL 




yRUlVR 








M v 













yRU2VR 








M c 













yRU3VR 








M T 



13 



M 



V 












UldiVl 




















VLD2VL 




















VLD3VL 


VrdiVr 








M D 











UrdzVr 








M s 











VRD3VR 








M B 



(60) 



It may be useful to comment on the diagonal forms of the matrices in Eq. fl57J|) . These simple 
block-diagonal forms follow from neglecting the off diagonal entries of the initial 6x6 quark 
mass matrices as in Eq. QB0D . Therefore, the forms of Eq.(^) become more complex in general 
if the off diagonal entries of the triangular matrices denoted by jjl and dr of the quark mass 
matrices are included. Because the tree level FCNC is already suppressed, in one loop level 
calculation, we can set Zij = 0. This is equivalent to neglecting the off-diagonal element of 
the triangular matrices yij in the one-loop calculation. Therefore, it is sufficient to keep only 
the diagonal element of y for the present purpose. In the approximation, there is a convenient 
parameterization of the 6x6 generalized KM. First the matrices shown in Eq.(^) can be 
diagonalized by the 2 by 2 block diagonal matrices, because we only need to diagonalize a 
singlet-doublet mixing in each flavor, 



Volu 

VoLD 



Clu 


Slu 


—Slu 


Clu 


Cld 


Sld 


—Sld 


Cld 



Vr 



ORU 



Vr 



ORD 



Cru 


Sru 


—Sru 


Cru 


Crd 


Srd 


— Srd 


Crd 



(61) 



where, 





( s Lu 








Slu = 




SLc 












SLt 




( Sr u 








Sru = 





SRc 







I o 





SRt 




( S L d 








Sld = 





SLs 







V o 





SLb 




f SRd 








Srd = 





SRs 







I o 





SRb 












Clu = 





CLc 






v 





CLt J 




/ C Ru 








Cru = 


to 


CRc 












CRt 









o > 


Cld = 





CLs 







v o 





CLb J 












Crd = 




CRs 







o 





CRb 



(62) 



with slu = sin 9lu and cl u = cos 9l u etc. and likewise for the down-quark sector. The flavor 
mixings of left handed quarks are given by the following 6x6 matrix: 



V 1 



Clu 
Slu 



—Slu 
Clu 



U L 




Cld Sld 
—Sld Cld 



CluU l Cld CluU l Sld 
SluU l Cld SluU l Sld 



(63) 
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where 3x3 part of the generalized KM matrix V L is no longer unitary, while U L is an unitary 
matrix. For the right-handed quarks, the CKM matrix is given as follows: 



yR 



CruU R Crb CrjjU R Sre> 

SruU r Crd SruU r Srd 



(64) 



We give an approximate parameterization of V L and V R . We first note that the mixings su 
and SRi are expressed in terms of the mass matrices elements: 



SRv 



VluiVl 
Mu ' 

VruiVr 



M, 



SLc ^ 



SRc 



u 



VLU2VL 

M c ' 

VRU2VR 

M c ' 



SLt 



sin 



c Rt = sin 



tan 1 



tan 1 



2\y LU3 \r) L M T 


\URU3 


Wr + M 2 - 


VLV3 


2 vl. 



2\Vru3\vr m t 


\llRU3 


Wr-m*- 


VLU3 


2 Vl 



(65) 



(66) 



There is a simple relation among the mixing angles, the light quark mass m qi and heavy quark 
mass rriQi which follows from zeros of the 3x3 ordinary light quark sector of the mass matrix, 

CRiC Li m qi = -SRiSuuiQi. (67) 

These relations are useful to estimate the box contributions in FCNC. We observe the following 
facts on the mixing angles between the singlet and doublet quarks. 

• The mixing angles are suppressed as su = 0{rrii/M WR ) and s Ri = 0(rrii/M WL ) for light 
five quarks (from up to bottom quarks )and their singlet partners. 

• s Lt is at the most O(f^) with M T = 0(rj R ). For M WR > 400 (GeV), s Lt may be less 
than 0.2. 

• SR t is not suppressed. In section 2, we argue that when Mt < Mwr, we can prevent 
the seesaw mechanism to act and the top quark mass is at the electroweak breaking 
scale without introducing non-perturbative Yukawa coupling. This corresponds to CR t = 
O(-p^) < 1 and s Rt is close to one. 

Because U L is 3 x 3 unitary matrix and we have direct measurements on and V c f, we can 
parameterize U L with the Wolfenstein parameterization. About the right-handed KM, we are 
guided by theoretical assumption. We assume that the Yukawa coupling between the doublet 
and the singlet quarks is left-right symmetric, = yR. This reduces to U L = U R . In the 
following section, we always assume this relation. It is instructive to write the generalized KM 
V L and V R explicitly in the approximate parameterization. We are interested in the charged 
currents between light down type quarks (d,s,b) and up type quarks. By neglecting the small 
mixings to the heavy quarks (U,D,S,B,C) ,4x3 part of the V L and V R are: 



( u L c L t L T L ) 



2 



1 

-A - iA 2 \ 5 r] 
c Lt A\ 3 (l -p-irj) 
\ s L ^A 3 (l - p-irj) 



A 



c L t{- 

SLt{- 



2 



1 - 

-A-iAX^X 2 
-A-iAX 4 r])X 2 



AX 3 (p - 
AX 2 

CLt 
SLt 



irj) ^ 



I d L \ 

sl 
\b L J 



(68) 
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( u R c R t R T R 



1-^ 

x 2 



A 



1 - 



A 2 



-A - iA 2 \ 5 rj 
c m AX\l-p-iri) c m (-AX 2 -iAX 4 ri) 
\ s Rt AX 3 {l -p-irf) s Rt {-AX 2 - iAX 4 r]) 



AX 2 

c Rt 
sat 



J 





( d R \ 




SR 




[b R J 



(69) 



Note that V L and V R is not left-right symmetric. In the limt of M t r = 0, they are given by, 



V 



1 - 



2 



-A - tA 2 X b r] 



X 

2 



AA 3 (1 -p-irj) -AX 2 - iA 2 X A r, 




V 



R 



1 



2 

|2\5„ 



-A - iA 2 X b r] 





A 

x 2 



AX 3 (p - irj) \ 
AX 2 
1 

J 

AX 3 (p - irj) \ 
AX 2 


1 / 



(70) 



(71) 



V AA 3 (1- p-iri) -AX 2 -iAX A rj 

6 One-Loop Level FCNC 

In the present model, there are many particles which may contribute to the FCNC in the 
one-loop level. For example, there are new contributions to the AF = 2 transition from 
box diagrams involving ordinary quarks and heavy isosinglet quark intermediate states. They 
can contribute to the Feynman diagrams in which WlWl, WlW r and W R W R are exchanged. 
W R W R exchanged diagrams are suppressed and we ignore their effect. The major contribution 
from WlWl and W r Wl exchanged diagrams is discussed below. 

In the WlWl exchanged diagrams, the lightest singlet up type quark T can contribute. We de- 
note Mi2(SM) as the contribution to the off-diagonal matrix element of neutral meson systems, 
i.e., KK and BB in the SM. Mi 2 (LL) is the WlWl exchanged box diagram in the present 
model in which T quark loop is taken account of. The deviation from the SM due to T quark 
loop can be written with Inami-Lim functions, 



M l2 {LL) -M 12 {SM) 



M 12 (SM) 



Af 



(Af L ) 



Mi 



,S(xt, Xt) 
' S(x t ) 
(k, I) = (d, s 

x T 



V4At + 



■ S(x T ) 



S(x t ) 



Ds 



Lt ) 5 



rriT 



{d,b),{s,b), 

\2 



WL 



'Mi 



(72) 



WL 



where (k,l) = (d, s),(d, 6), and (s, b) correspond to KK , BjBd and B S B S mixings respectively. 
S(x) and S(x, y) are the Inami-Lim functions, su is the mixing angle of the singlet quark T 
defined in Eq.(|6~5D. The mixing angle is approximately given by, 



SLt = O 



M WL M 7 
M WR 



(73) 
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We find that for My = O(Mwl), the mixing angle is suppressed as su = 0( M Y L )- For Mr = 

M WR 

0(M WR ), we obtain the larger mixing angle, i.e., s Lt = O(-p^). When M WR ~ 400(GeV), 
for the former case, the mixing angle is 0.04 and for the latter case, it is 0.2. Therefore, for 
the smaller M T = 0(M WL ), the T quark does not contribute to WlWl exchanged diagrams 
at all. For the larger mixing angle case, the deviation from the SM can be as large as 30% 
for Mwr — 400GeV. In the subsequent analysis, we extensively study the case for the small 
mixing angle. Then WlWl box diagram contribution is not changed from that of the SM, i.e., 
Mi 2 {LL) ~ M 12 (5M). For the latter case, we need to do more careful analysis and this will be 
done in the future publication. 

When the T quark contribution to WlWl exchanged diagrams is small, the most important 
contribution comes from WrWl exchanged diagrams. |TI| In the WrWl exchanged diagrams 
(Fig.l), we must also take account of both isosinglet and isodoublet quarks as intermediate 
states. 



6.1 K° — K° mixing 

The K° —TC° mixing is one of the best processes to test the FCNC in our model. We reexamine 



the famous enhancement factor in Wl — Wr box diagrams of LR models |19] taking account 
of both isosinglet and isodoublet quarks as intermediate states. The effective Hamiltonian for 
Wl — Wr exchanged diagram is, 

C 2 

Hes(LR) = -^M^ L A^ R Af L 2/3^E^F(x y , x s , (3)(dLs)(dRs) , (74) 



where 
and 



Af = V£V* , Af = V q f d V q L ss , q^ s = («, c, f , U, C, T) , (75) 



/?= — —) , F(x 7 ,x 5 ,/3) = (A + (3x 7 x s )Ii(x 7 ,x s ,(3) - (1 + f3)I 2 (x 7 ,x s , /3) , (76) 



M Wr 



X ^ X ^WJ ^ =1 ~ 3 >- x ^ x ^ = WJ «^ 4 ~ 6) - (77) 

where i runs from 1 to 3 and V L I R denotes the 6x6 mixing matrices. The loop functions I\ 
and I 2 are found to be the same ones given by Ecker and Grimus [[20H : 



Xry In cc^i / \ $ ■ — ■) ■ 

I 1 {x^x s ,P) = 7 - v , w^j + (x 7 <-> x 5 - — — — — , 78 

(1 — X 7 )(X 7 — X5)(l — /5X 7 ) (1 — PX 7 )(1 — pXs){l — p) 



111 x In /9 

J 2 (x 7 , x 5 , /?) = 7 7 r + (x 7 «-> xa) - — — — — — . (79) 

(1 — x 7 ){l — px 7 )(x 7 — xs) (1 — p)(l — px 7 ){l — pxs) 

The contribution to the off diagonal matrix element M^ is 
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M«(LR) = ^^M^ L A^ R Af L 2P y /x^xsF(x J , xg, ft) (K°\(dLs)(dRs)\K°) (80) 
The Ki and K$ mass difference, Amx and the CP violating parameter ex are given: 

Am^ = 2ReM^, = exp(z-)( 12 + ft , £o = p-r^. (81) 

The Li? contribution to Amjf is given by 

Am K (LR) = 0M^/>^ ]T Re(A^ L )2/5 v ^-r7F(x 7 ,x (5 ,/5) , (82) 

7,5=1 

where 

= ^Kf 2 K m K , (83) 

with 

4(m s + m d ) 2 8 V ; 

We estimate the matrix element using the vacuum-insertion approximation. In the following 
estimate of the LR contribution, we neglect the effect of the QCD corrections. 

The sum over 7, 6 = 1 ~ 6 can be reorganized in a better way. By using eq. (|6TD , we obtain 



c Lqi c Rqiy /x q ~ = -s Lqv s Rq ^X Qi , (i = 1, 2, 3) (85) 
which lead to a useful formula 

Am K (LR) = ^M^J 2 K m K n2px 
3 

Re(\f R \f L )c Lqi c Rqz c Lq] c Rqj ^x^-F(x qi ,x qj ,X Qt ,X Qj ,P) , (86) 

where 

F(x qi , x qj , X Qi , X Q . , /5) = F(x qi , x ?j , /?) - F(X Qi , z Q , , (3) - F(x qi , X Qj , (3) + F(X Qi , X Qj , /?) , 

(87) 

which is a gauge invariant quantity [21]. Af R is defined by 

\f R = U q %U q R s , (1 = 1,2,3). (88) 

Before closing this sub-section, we show the SM contribution Mf^(LL), which is used in the 
following calculations: 

M*(LL) = ^F 2 LVm K M^[(Af ) 2 r ?1 5(x c ) + (X^ L ) 2 m S(x t ) + 2X^X^ L r ]3 S(x c ,x t )}, 

\ LL = U%pZ., (i = l,2,3), (89) 

where S'(x) and S(x,y) are the Inami-Lim functions, and 771 = 1.38, 7/2 = 0.57 and 773 = 0.47 
are used. The parameter E>k is taken to be 0.75 ± 0.15 in the following calculations. 
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6.2 Numerical results in LR symmetric limit 

In the limit of m u = 0, the up flavors contribution, i.e., the (i=u and/or j=u) contribution to 
AmxiLR) vanishes. The loop integrated functions F(xi,Xj,Xi,Xj, (3) for charm flavors, top 
flavors and the mixed flavors intermediate states are given approximately as: 



F(x c , x c , X c , X c , 0) ~ 5 + 41nx c - \ny c + In (3 

A-x t x 2 t -2x t + A 
I- x t {I- x t 



F(x t ,x t ,X T ,X T ,P) ~ - + ^ A^-ln^ + ln/? 



4 _ T 

F(x c ,x t ,X c ,X T ,f3) ~ ^-Jlnxi + ln/3, (90) 

1 Xt 

where jji = = rriq./M^. The contribution of each intermediate state also depends on the 
mixings: 

Xf R Xf L = U^U R s U R * d U^ s . (91) 

In the left-right symmetric limit of the Yukawa couplings, i.e. Vl{u,d) = Dr{u,d) hi eq. (|60|) , 
the relations U^ r , = C/fL and = U R hold. As we discussed in the section 5, we can 
apply the Wolfenstein parameterization on U L and U R . Taking into consideration these mixing 
matrices, we can estimate contributions of charm flavors, top flavors and mixed flavors of 
charm and top flavors intermediate states. The charm flavor contribution includes (c,c), (c,C), 
(C,c) and (C,C) as intermediate states. The top flavor contribution includes (t,t), (t,T), (T,t) 
and (T,T). The mixed flavor contribution comes from (c,t),(C,t),(c,T), and (C,T) intermediate 
states. The relative contributions of ordinary quarks and singlet quarks are found by estimating 
F(x qi ,x qj , Xq^Xqj , (3). Using the Wolfenstein parameterization, the mass difference Am^(L_R) 
is given by, 

AmK{LR) = 2V2C e K[3{\ 2 x c F(xc,x c ,X c ,X c ,P) 

+ A 4 A 10 {(1 - p) 2 - r, 2 }x t c 2 Rt F{x t , x t , X T , X T , (3) 
+ 2A 2 \\l - p)^x~x~ t c Rt F(x c , x t , X c , X T , (3)}, 

Ce = = 3, 78 x 10 4. (92) 

6\'2TT 2 Am K 

x c F(x c , x c , X c , X c , (3) ~ -0.015, x t F{x t , x t , X T , X T , (3) ~ -14 and ^/xjT t F{x c , x t , X c , X T , (3) ~ 
0.018 for Myyn = 2(TeV). So, it is found that the charm flavor contribution is the most im- 
portant in AmxiLR) .We also note that the Wl — W R box graph contribution to Am^- is 
negative relative to the Wl — Wl one. Among the charm flavor contributions, a quarter of the 
ArriK(LR) comes from two isosinglet quark intermediate states, i.e., (C,C). This can be seen 
that the singlet C quark intermediate states contribution in F(x c , x c , Xc, Xc, (3) is approxi- 
mately given by 1 — \ny c + \n(3. In the numerical calculation, we must give the mass of the 
isosinglet quark C. It is determined as follows. The mass of charm quark is given by the seesaw 
formulae, m c (W m c)?/i?[/2?/Li/2?7L- In the left-right symmetric limit, i.e., y RU2 = y RU 2 = Vm, 
it becomes, m c = 2 ( ^ ) 2 Mw ^wl _ \yith the assumption of the strength of Yukawa coupling, 
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UU2 = 0(g), we get mc = 0(100 x Mwr)- In our numerical calculation of ArriK and ex, we 
set m c = 100M Wr . The Wl — Wl box graph contribution is well known and is given by 



= V2C e B K \> XcVl . (93) 

Am K 

The W L — W R box graph contribution to Am K is negative relative to the W L — W L one. 
Therefore the experimental value of K L and K s mass difference can be fitted only if there is 
sizable constructive contribution to Am K due to the long distance effects or if W R is sufficiently 
large and ArriK is saturated by the standard model like contribution. Including possible long 
distance effects, ArriK consists of three parts. 

Am K = Am K (LL) + Am K (LR) + Am Long , (94) 

where Am Long is the long distance contribution. First, we show the ratio ^j^ff j versus M Wr 
in Fig.2, where m s = 120MeV or 200MeV, and B K = 0.75 are used. In orderW the W L - W R 
contribution to be smaller than the W L — W L one, M WR > 1.3TeV should be satisfied. In Fig. 3, 
we show the AmLong which is needed to explain the experimental value Amx in the present 
model. If we allow the sizable Am^ ong = O(ArriK), Mwr is allowed to be ITeV. Therefore 
we conclude the lower bound for Mwr obtained from Amx is O(lTeV), though there is large 
theoretical uncertainty on the lower bound on Mwr from Kl and Ks mass difference. 



6.3 CP violation 

Let us discuss the CP violation of the K° — K° system. By studying €k in the present model, 
we can obtain alternative constraint on two important parameters M W r and M T in the model. 
As we discussed in the previous sub-section, the lower bound on M W r is obtained from Am K . 
However it is shown that and Ks mass difference is not sensitive to M? because the top 
flavors contribution is tiny in the real part of M 12 . However, in the imaginary part, the top 
flavors contribution becomes important as Mt is larger. Neglecting tree level FCNC, there is 
one CP violating phase for the left-handed mixing matrix U L and one for the right-handed 
mixing matrix U R . In the left-right symmetric limit, U L = U R holds and we have only one 
CP violating phase. The CP violating effect of the Wl — W R exchange in the K° — K mixing 
is proportional to the CP violating phase. The imaginary part of M^, which comes from 
Wl — W R exchanged diagrams is given as follows: 

lmM«(LR) = ^L MwL f 2 K m K K2p x 

3 

]T lm{Xf R Xf L )c Lqi c Rqi CL q3 c Rqj ^x^-F{x qz ,x qj ,X Qi ,X Q] ,(3) (95) 

where 

lm(\L R \ RL ) = -2A 2 X 6 V , Im(\^ R \ RL ) = 2A 4 X 10 r](l - p), Im(A^Af L ) = 2A 2 X 6 r]. (96) 
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€k in the present model can be written as, 



(-LL 



e LR = C e K(3rj{-2A'X b x c F(x c , x c , X c , X c , (3) + 2A 4 A iU (l - p)c l m x t F{x t , x t , X T , X T , (3) 



7I\ 



+ 2A A y/x c x t c Rt F(x Cl x t , X c , X T , (3)} exp(i-), 



1 



B K C e ri{-^ 2 \%S(x c , x c ) + 2A 4 A lu r/ 2 (l - p)S(x t , x t 



4\10„ 



7T , 



+ 2A 2 \%S(x c ,x t )}ex P (i-) 



(97) 



For numerical calculation, we must know how Cm and tjit depend on the parameters Mr and 
Mwr. In the present model, the masses of the top quark and its singlet partner are given by, 



rriT 



Mi 



1 



V 9' 



w L - 



I V2M W 



V2M } 




M> 



T 



\V2M ] 



(98) 



(99) 



w R , 



(c.f. Eq.(H)). In the case of the left-right symmetric limit i/l = Ur = y, the Yukawa coupling y 
is given in terms of the top quark mass: 



y =9 



2 m t 
4M^ 



\ 



w R m t 



(100) 



Then the singlet quark mass is given as 



rriT 



Mw R 1 , , 



Mi 



\ 



1 + 4 



M^ £ M | 
M^mr 



(101) 



In the case of M Wr ^> M T , we get m T ~ miM^/M^. Thus, the physical mass of the 
singlet T quark is determined if My/ R and Mt are fixed. In Fig.4, we show the magnitude of 
Yukawa coupling in the region of Mt < My/ R - The Yukawa coupling increases as Mt is larger. 
Increasing Mt further may invalidate the perturbative calculation. In Fig.5, crj is shown versus 
M t /M Wr . As far as M T -C M Wr , the singlet-doublet mixing c^ t is approximately given by 
CM ~ M t /M Wr . Coming back to €k, the top flavor contribution is suppressed by the factor. 
Such suppression mechanism is distinctive in the present left-right model, cm increases as Mt 
is larger. So the top flavor contribution is more important as Mt increases. 

There are two distinctive cases on the contribution to e parameters. One corresponds to the 
case that is positive and the other corresponds to is negative. The ratio is shown as 
a function of Mt in Fig.6, where Mw R = 0.5, 1.5, and 3 TeV are taken. The ratio does not 
depend on rj and we set p = 0. As seen in this figure, the LR contribution is negative in 
M T > 60GeV, but positive in M T < 60GeV relative to LL contribution. (See Eq.@7p.) If 
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£lr 7^ and we would use the standard model expression to fit ex, i.e., with e# = ^ll(p,v)i 
there would be disagreement between the allowed region of (p, 77) determined by B physics 
data only and that obtained by ex- Specifically, here we mean \V u b\, \V c b\ from semileptonic B 
decays and \V t d\ from BB mixing as B physics data. When we obtain the constraint on \V t d\ 
from the BB mixing, we can use the SM expression for M 12 because the deviation from the 
SM in WjJWl exchanged diagrams are negligible for the small mixing angle sit = 0(( ^^ ) 2 ) 
and WrWl exchanged diagrams do not contribute to BB mixing as we show later. After all, 
for the small mixing angle case, the B physics data is not affected by the presence of the new 
physics. Then the mismatch between the constraint of ex and that of the B physics data will 
be manifested in the following ways. 

(1) The (p, 77) determined by the B physics data is above the line obtained from the constraint 
€k = eLL(p,r]). This case corresponds to €lr/€ll < 0. 

(2) The (p, 77) determined by the B physics data is below the line obtained from the constraint 
€k = £ll(p, if). This case corresponds to e LR /en > 0. 

Qualitatively, this mismatch is understood as follows. If the CP violating parameter 77 deter- 
mined by the B physics data is larger than the one determined by the standard model constraint 
€k = we need some negative contribution to ex relative to en- If the 77 determined by 
the B physics data is smaller, then the en is not enough to explain ex and we need positive 
contribution to ex- The present data of B^Bd mixing and \ V u b\ as well as the bag parameters 
Bx and Bb are not sufficiently accurate to discriminate the one case from the others. Of 
course, this mismatch may disappear by using the full expression for ex in the present model 
with the parameters M T and M WR suitably adjusted. In Fig. 7, we show the allowed region for 
{Mwr, j\f~) using the full expression for ex, 

ex = e LL (p, 77) + e LR (p, 77, M WR , M T ). (102) 

Taking into account of the present allowed region for (p, 77) obtained from BdBd mixing, |V^| 
\Vyjb\ and ex, the allowed region for (M WR , ^ ) is plotted. Depending on the two cases, 
the allowed regions of (M WR , ^ ) are quite different from each others. The lower curve 
corresponds to the case (2), i.e., ^ > 0. The upper bound on Mwr is obtained in this 
case. The upper curve corresponds to the case (1), i.e., ^ < 0. Since 77 is larger than 
one which would be obtained by the SM fit, the negative contribution is needed. This is 
achieved by taking larger. To maintain this effect while taking M WR larger, M t r must 
increase. We use the following experimental inputs to determine the allowed range of (p, if) : 
|yCKM| = 3 95 ± 0.0017, m t = 175.6 ± 5.5GeV, \V£ KM /V$ KU \ = 0.08 ± 0.016 and |V^ KM | = 
0.0084 ± 0.0018 13 @. We also use B K = 0.75. The above constraints are shown in (p, 77) 
plane. (See Fig. 9.) In order to find the Bx dependence of our result, we also show the allowed 
region by taking Bx = 0.6 ~ 0.9 with p = and 77 = 0.35 in Fig. 8. We can also discuss the 
Wl — W R exchange effect on B® — B d mixing. The contribution to Am B o is given as 

Am B o(LR) = p- 2 M w J 2 Bd m B K2f3 x 
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3 

^ \\^ R \f L \c Lqi c Rqi c Lq] c Rq]y /x^F(x qi1 x qj1 X Qi1 X Qj1 f3) , (103) 

where k ~ 3/4 and 

A^Af = U^U^U^U^ . (104) 

In the SM, the top quark intermediate state dominates the B d — B° d mixing. The Wl — Wr 
exchange contribution is also dominated by the top flavor intermediate state, however, this 
contribution is suppressed by CR t , that is of order 10~ 4 compared with the SM contribution. 
Thus, for B d — B d mixing, the LR effect is negligible. 



7 Rare Decays of K Mesons 

Experiments in the K meson system have entered a new period with the observation of the rare 
process K + — > 7r + vV, and the dedicated search for Kl — > 7r°z/z7. Recently, the signature of the 



decay K + — > n + i/V has been observed by E787 Collaboration [15] and the reported branching 
ratio is 4.2lg; 5 x 10~ 10 , which is consistent with the value predicted by the SM. Additional 
(and improved) data are expected in the near future. In view of this situation, a detailed 
study of the rare K meson decays is necessary. The decay K L — > n^vV is one of the most 
promising processes, since it is a CP violating mode in the SM. This mode is theoretically 



clean to extract the CKM parameter 77 |18[ . We investigate rare K meson decays in the present 
model introducing right handed neutrinos. However, in the model, the neutrino masses are 
zero in the tree level and lepton flavor is well conserved ( see analyses within other models 
|24||). Scalar and tensor operators appear due to the LR box diagrams, in which both the left 
and right handed gauge bosons, Wl and Wr, are exchanged. The scalar operators have an 
enhancement factor MK/m s in the matrix element < Tr\sd\K >. Thus, the scalar operator may 
make a large contribution to the rare K meson decays, K + — > tt + vV and K^^ ~ > ^vv. An 
important point is that the CP property of the scalar interaction is different from the V-A 
interaction (sd)y _ a^Pi v i)v - a in the SM. The decay K L — > n i/V through the scalar operator is 
not a CP violating one, so we have a non-zero branching ratio B{K L — > it°viVi) even in the CP 
conserved limit (rj 0). Thus, it is important to estimate the size of the effect of the scalar 
operator on the pion energy spectrum. 

7.1 Rare decays by scalar operator 

The rare K ttuV decays are loop-induced FCNC processes in the SM, that, being dominated 
by short-distance physics, are theoretically very clean modes ||18|| . The matrix elements involved 
in these decays are related to the experimentally well known decay K + — > 7r°e + z/ using isospin 
symmetry, and corrections to this relation have been studied . We start our analysis of such 



decays with the following effective Lagrangian, which is produced from the LR box diagram in 
addition to the SM contributions: 

= J2 [ C l SM (s L ^d L )(V wl ^ Ld ) 

V 1 l=e,fi,r 
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+S l LR {s L d R )(uL,i^R,i) + S l RL (s R d L )(v Rt iv Lt i) 

+T l LR (s L a^d R )(v Ljl a lxu u R>l ) + T RL {-s R a^d L ){v R ^^{) + h.c] , (105) 



where k = a/(27rsin <dw)- The first term is the SM contribution [18, EBf, and from the second 



to the fifth are the new contributions arising in this model. The scalar and tensor operators 
generally appear from box diagrams when one considers a model which contains the right- 
handed charged gauge boson W R . As a concrete analysis, in section 7.3, we investigate the pion 
energy spectrum for LR symmetric scenario, S l LR = S l RL and T l LR = T l RL , which correspond to 
the limit cos6> Ls = cos9 Rd = 1 with U R = U L . 

There are also penguin diagram contributions to the process we are interested in. However, 
only the box diagrams produce the scalar operator (sd)s(yv)s m the effective Lagrangian, thus, 
we do not consider the contributions from penguin diagrams in this paper. 

First we show the decay amplitudes for the neutral K meson states and K$ Q. The decay 
amplitudes A(K L S — > vr°z7z/) are: 

A(K L , S -> 7r°F M ) = ((pC l SM T qC l s * M ) < Vfd > (P, 7m (1 - 

+ (P(S 1 LR + S l RL ) ± q(S l LR + S l RL y) <sd> {v lVl ) 

+ (p(s 1 lr - s l RL ) =f q(s l LR - s l RL y) <sd> {vasvi) 

+4 (p(T l LR + T l RL ) T q(T l LR + T l RL )*) < sa^d > faa^) 

+4 (p(T l LR - T l RL ) ± q(T l LR - T l RL y) < sa^d > (u^^)) (106) 

The CP conserved limit corresponds to p = q, with all coefficients C l SM , S l LR , S l RL , T l LR and 
T l RL real. In this limit, the decay amplitude A(Kl — > tt°I7 l z/ l ) through the V-A interaction is 
zero, while A(K$ — > tt^l^l) is nonzero, and the decays through the scalar operators A(Kl^ 
^VhVlj^VlVr) remain non-zero generally. In the LR symmetric parameterization, Kg decay 
through the scalar operators is the CP violating mode, while CP is conserved for Kl decay. 
Thus decays of neutral K meson in the LR symmetric parameterization are summarized as 
follows: 

{sd) V -A{viVi)v-A C/P 
•i^decay { ^sd)s{pii l i)s =>• CP Conserving , 

{sd) T {vm) T => c/p 

(sd) V -A(viVi)v-A CP Conserving , 
»K S decay { {sd) s {Vivi) s => C/P 

(sd)T(ViVi)T CP Conserving . 

Experimentally we do not observe the neutrinos, and the pion energy spectrum is obtained 
by summing these contributions which have different CP properties with each other. The 
Kl decay through the V-A operator is suppressed due to CP symmetry, but decays through 



1 We use a conventional phase choice where |-K"z,,s >=p|i^° > ±q\K° >, CP\K° >= — \K >. 
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the scalar operators are CP conserving ones, furthermore their matrix elements are enhanced, 

< Tr°\sd\K >~ MjLf± ( where < (sd)v-A >= f+P+ J rf-P--,P± — Pk^P-k ), as we can see from 
the equation of motion in the next subsection. Thus, the contribution of the scalar interaction 
to the decay amplitude A(K L — > ifivv) is sizable and dominates in the CP conserving limit. 
The decay amplitude for the charged K meson, A(K + — > 7r + z/z7), is obtained in the same way: 

G 

A(K + -> ti+Vm) = ( °SM < srfd > {V lllx (l - 75)^) 

+ (S LR + S RL ) <sd> (Vm) 

+ (S L r ~ S RL ) <sd> {Va^vi) 

+2 (Tlr + T RL ) < sa^d > {V l( j^n) 

+2 (T LR - T RL ) < sa^d > (!7,<V75i/,)) , (107) 

where < O >=< ti + \0\K + >. 



7.2 The matrix elements and the coefficient functions in the LR 
model 

In this section, we explain our estimations of the matrix elements and show explicit forms of 
the coefficient functions in the LR model. The matrix elements in the SM, < (sd)v-A >, can 
be related to the matrix elements of the experimentally well known leading decay K + — > 7r°e + z/ 
using isospin symmetry. So the form factors f^ + - >7r+uu an d f^°^^ vv are written in terms of 
j?K ^n°e+v_ rjj^ ma ^ r j x e i emen t f the scalar operator < (sd) s > is also related to these by the 
equation of motion: 



< n\sd\K > 



p_- < n\(sd)v-A\K > 
m d - m s 

= U^^ + f-^—, (108) 
m d -m s m d - m s 



where f± are the form factors of the corresponding matrix element < -K\{sd) V -A\K >. We 
estimate the matrix element of the tensor operator using the NJL model. First we note that 
the matrix element of the operator is the second order of the typical momentum. Therefore 
in the sense of the chiral perturbation , it is enough to estimate the form factor fa at zero 
momentum transfer, which is defined by < sa^d >= ■ (p+p 1 ! —P-P'+j- the SU(3) breaking 
effect may be also neglected. Calculating the Feynman diagram shown in Fig. 10, the form 
factor is given by 

where m q is a constituent quark mass and f n is the pion decay constant. As a rough order esti- 
mate, we assume m q = 200 MeV, f w = 120 MeV. Then we obtain , f T = — 0.37i. Substituting 
the value in the pion energy spectrum, we can see that it is safely neglected. 
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Now we show the explicit form of the coefficient functions S l LR and T l LR in the model. We 
calculate box diagrams, in which the left handed Wl boson and the right handed gauge boson 
W R are exchanged, as seen in Fig. 7. There are corresponding charged Higgs diagrams due to 
the gauge invariance. The internal upper fermion lines correspond to the ordinary and the 
singlet quarks, the lower ones correspond to the SM and the singlet leptons. 

The coefficients in the effective Lagrangian are: 

S 1 lr = CL.au E (U^U* d )c Lq c Rq x{3^x-mF(x q ,y h X Q ,Y L ,{3), (110) 

q=u,c,t 

Tin = -asCRd E (UqsUqd)aqCRq x (3^/XqljiF' (x q , yi, Xq, Y l , (3) . (Ill) 

q=u,c,t 

where 0L/Rq is the mixing angle between a singlet left /right handed quark and the corresponding 
doublet quark as defined in Eq.(|^), and UJ are 3x3 CKM matrix elements. In addition to 
the parameters x q ,Xq and (3 of Eq. (|77|) , new dimensionless parameters Vi,Yl are defined by: 

2 2 

mi m T 
Vi = Y L = jjt, (112) 

m w L m w L 

where mi are mass eigenvalues of the ordinary leptons e, \i and r, while rriL are the corresponding 
additional heavy lepton masses. 

The function F ' is defined by replacing F(x,y,f3) with F '(x,y,/3) in Eq.(|87f) : 



F'{x,y,P) = -h(x,y,/3)-^I 2 (x,y,/3). (113) 

The other coefficients, S l RL and T RL , can be obtained by interchanging the indices L and R. In 
the limit when cl, s = c,R,d = 1, S l LR = S l RL and T l LR = T RL . Thus, we simply write them as S l 
or T l and take the coefficients as LR symmetric. The coefficients S (I = e, fx) for electron and 
muon are negligibly small due to the smallness of y e and y^, and only S T contributes to the 
process significantly. The coefficient function S T evaluated with M R = 500 GeV is: 

S T ^ -6.10 x 10" 7 A M - 8.76 x 10~ 5 A C - 1.46 x lQ- 3 c Rt \ t , (114) 

where X u ~ A (l — 4r) , A c ~ — A (l — ^) an d A t ~ — A 2 X 5 (1 — p — irj) in the Wolfenstein pa- 
rameterization. The coefficient for the up quark is negligible. For the top quark, the coefficient 
function is enhanced by its heavy mass, but suppressed by the CKM factor At, as compared to 
A c from the charm quark. Thus, the contribution from the top quark in the model is small in 
the case of LR symmetric case. It is important that all the coefficients of the LR model are 
suppressed by the factor /3, so the V-A interaction of the SM dominates as M R becomes large. 



7.3 Pion energy spectrum 

We present the pion energy spectrum obtained using the coefficients of Eq. (|110|) and Eq.( |lll| ). 
In the process Kl — > 7r°z/77, the contributions from scalar interactions are controlled by (3 and 
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the enhancement factor M^/m s in the matrix element ((sd)). For SM, the decay amplitude of 
Kl — > 7r°z/z7 is proportional to rj, thus the contributions of scalar interactions relative to that of 
SM for the pion energy spectrum is O (J^^J • In the process K + — > n + vV, contributions from 
the scalar and tensor interactions are tiny compared to the SM one. The pion energy spectrum 
for the decay K + — > n + vV is: 



dB (K- 



5 N ■ k + ^xI-A5 2 



x 



(xl-45 2 ) (\C l SM \ 2 + l(kvli\Tf)+3i 



K 



m d - m t 



(l - 5 2 + £ if \S l \ 



k+ = 
5 N = 



2n 2 sm 9 W 
\f+(K+ ^it°e + u)\ 2 

Jdx„\f+(K+ -> 7r°e+z/)| 2 (xl - 45 2 ) 3/2 ' 



(115) 
(116) 
(117) 



where normalized pion energy defined by x^ = 2E 7T /Mk, and 5 = m^/Mx, £ = 

/r^/zr* . ^ = \fr + *v/r"°i and * = (! + 52 - 

The process fTf, — ► vr°z/z/ is a more sensitive probe of the scalar interactions, as discussed in 
the previous section. The energy spectrum is given by: 

dB [K L/S 



X 



= 5 N ■ K L ^jx\- A5 2 
xl - A5 2 ) ( \pC l SM =f g^ M | 2 + co 2 t \pT l T q T l *\ 

(l - 5 2 + £ if \pS l ± qS l * 



+3t 



M, 



m d - m t 



(118) 



where cuq 



K°7T° I 



The first term is contributions from V-A 



and. L = 

and tensor operators, SL,R<J^ v dL,R, the second is the contribution from the scalar operator 
which has an enhancement factor M\j{md — m s ) 2 in the matrix element. The contribution 
from the tensor operator is suppressed by the kinematical factor (x 2 — A5 2 ) in the low pion 
energy region, x n ~ 25, whereas at large pion energies, x n ~ 1 + 5 2 , it is suppressed by the 
factor t. Furthermore, decay through the tensors is CP violating and has no enhancement 
factor. Hence, the contribution from the tensor operator is negligible compared to the scalar 
operator and the SM contributions. The form factors are related to those of the experimentally 
well known decay mode K + — > n°e + v. In Fig. 12 we show the pion energy spectrum, dB J^^ , 
multiplied by a factor 10 10 . To see the dependence of the branching ratio on rj, we plotted three 
curves, which correspond to the cases rj = 0.25, 0.3, 0.35. The solid lines are the pion energy 
spectra in the LR model, the doted lines are the corresponding SM prediction, where we have 
taken M^ s = 500(GeV) and p = 0.25. The LR contribution is large in the low energy region 
x n ~ 25, while in the high energy region the SM contribution dominates. In Fig. 13, we study 
the dependence of the pion energy spectrum on M WR . Notice that as the coefficient functions 
of the effective Lagrangian are proportional to (3, dB/dx w is proportional to (M Wl /M Wr ) 4 . The 
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effect is negligible for M WR > 1 (TeV) and the energy spectrum reduces to the prediction of the 
SM. 

The effect of the new physics can be seen in the pion energy spectrum of the Kl — > ifivipi 
decay for Mw R < 1 TeV. On the other hand, the analyses of K° — K mixing have given a 
constraint of M Wr > 1.6 TeV, which may discourage the search for new physics in this decay 
mode. However, it is important to comment on the constraint Mw R > 1-6 TeV. This bound 
has been obtained from the box diagram of the Wl — Wr intermediate state. We assumed, as 
is usual, that the short distance effect dominates K° — K mixing. However, if long distance 
physics dominates the mixing, the constraint of Mw R > 1.6 TeV is no longer valid. 

8 Conclusion 

We present the formalism and systematic analyses of the seesaw model for quark masses. 
A framework allowing the top quark mass to be of the order of the electroweak symmetry 
breaking scale is explained. The derivation of the quark mass formulae is presented in detail 
with flavour mixing included. There, it is shown that expanding simply by the inverse power 
of the singlet quark mass matrix fails and we propose an alternative expansion to overcome 
the problem. Furthermore, we find that a quark basis in which the singlet and doublet mixing 
Yukawa coupling is a triangular matrix is appropriate for finding the mass base. Starting in 
such a basis, by neglecting the flavor off-diagonal Yukawa couplings, we can reproduce the 
quark mass formulae which are obtained as the solutions of the eigenvalue equation. Also, we 
give the theoretical formulae for the tree level FCNC. The tree level FCNC in the model are 
naturally suppressed as (quark masses) 2 divided by an (SU(2) breaking scale) 2 . The effect on 
rare K and B decays is far below both the present experimental bound and the prediction 
of the standard model. As for FCNC beyond the tree level, the one loop effect involving the 
right-handed gauge boson exchange is discussed for K° — K mixing and for K — > ixvv. Arnx 
and ex give constraints on the parameters of the model Mt and Mwr- We reanalyse the 
Beall, Bander and Soni bound for Mwr in the present model and show the lower bound of 
Mwr is about O(lTeV). The constraint of ex is interesting and we show the allowed region in 
(Mwr, ^ t )■ If the allowed region of (p, 77) is tightened by the data of the B factory and/or by 
the improvement of the lattice computation of the hadronic matrix element, the allowed region 
will be much more specified than the present one. Alternatively the effect of the new physics 
can be seen the K L — > iPvipi decay in the case of M W r < l(TeV). There is a scalar operator 
in the effective Lagrangian, which come from LR box diagrams. For the decay Kl — > 7r°uV, 
there is a significant contribution from the scalar operator, especially in the low energy region 
of the pion energy spectrum, which, for the values Mw R = 500GeV and p, n = 0.25, amounts 
to an enhancement of about 30% to the total branching ratio. Thus, measuring the decay 
K L — > 7r°z/z7 precisely may be important to probe the effect from new physics. Other aspects 
of the present model, such as the constraints from the precision measurements and the other 
B decays, the Higgs sector, and the neutrino mass and mixings will be discussed elsewhere. 
The further improvement of the QCD corrections to our computation is also needed for serious 
comparison with the experiments. 
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(a) (b) (c) (d) 

s ds ds ds d 




Figure 1: Box diagrams for K° — K mixing involving singlet and doublet quark intermediate 
states. 
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Figure 2: ^|^| as function of M WR . The solid line corresponds to m s = 200(MeV) and 
the dashed line corresponds to m s = 120(MeV). We use Bk = 0.75. 
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Figure 3: The long distance contribution Am Long which reproduces the Kl and K$ mass 
difference in the present model. The solid line corresponds to m s = 200(MeV) and the dashed 
line corresponds to m s = 120(MeV). The unit is the experimental value of Am^. We use 
B K = 0.75. 
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Figure 5: The singlet and doublet mixing angle cos 9m as function of 
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Figure 6: ^ as function of M T . The solid line corresponds to M WR = 3 (TeV), the dashed line 
corresponds to Mwr = 1.5(TeV) and the dotted-dashed line corresponds to Mwr = 0.5(TeV). 
p is set to be zero. m s = 160 (MeV) is used. 
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Figure 7: The allowed region for (Mwr, Mwr ^' ^ e region between the upper curve and the 
lower curve is allowed. We take account of the present constraint for (p, 77) obtained from B d B d 
mixing and \V uh \. B K = 0.75 and m s = 160(MeV) are used. The upper curve corresponds to 
(p, rj) = (—0.087,0.43) and the lower curve corresponds to (p, 77) = (0.25,0.14). 
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Figure 8: The allowed region for (M WR ,^^) with (p,rj) = (0,0.35). The lower curve 
corresponds to B K = 0.6 and the upper curve corresponds to B K = 0.9. The region between 
two curves is allowed. m s = 160(MeV) is used. 
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Figure 9: The allowed region in (p, 77) plane obtained from \V u b\, \V c b\ and BdBd mixing. To 
extract \Vtd\ from BdBd mixing, we use the SM expression for Xd- This is a good approximation 

for the small angle case, i.e., su = 0( ( ^ I WL ) 2 )• (See also the text below Eq.([72|).) 
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Figure 10: The Feynman diagram for the matrix element of the tensor operator in NJL model. 
® denotes an insertion of the tensor operator sa^d. 




Figure 11: Box diagrams which contribute to the effective Lagrangian for the process K — > i^vv. 
(a) is a contribution from Wl and Wr. (b) and (c) are gauge boson and unphysical Higgs 
contributions, (d) is a contribution from the unphysical Higgs Xl and Xr- 
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Figure 12: The pion energy spectrum dB[Ki — > -K°vv]/dx for various values for the CP violating 
parameter 77. We use m s = 100(MeV), M WR = 500(GeV) and p = 0.25. SM denotes the 
predictions of the standard model and LR denotes the predictions of the left-right model. 
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Figure 13: The pion energy spectrum dB[Ki, — > ir°vu]/dx for various values of Mw R - th s = 
100(MeV) is used. SM denotes the prediction of the standard model and LR denotes the 
predictions of the left-right model. 
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